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^ ; ABSTRACT 

Aims. We analyse the influence of the stochastic particle acceleration for the evolution of the electron spectrum. We assume that all investigated 
^ spectra are generated inside a spherical, homogeneous source and also analyse the synchrotron and inverse Compton emission generated by 
such an object. 

Methods. The stochastic acceleration is treated as the diffusion of the particle momentum and is described by the momentum-diffusion 
equation. We investigate the stationary and time dependent solutions of the equation for several different evolutionary scenarios. The scenarios 
are divided into two general classes. First, we analyse a few cases without injection or escape of the particles during the evolution. Then we 
investigate the scenarios where we assume continuous injection and simultaneous escape of the particles. 

Results. In the case of no injection and escape the acceleration process, competing with the radiative cooling, only modifies the initial particle 
^ I spectrum. The competition leads to a thermal or quasi-thermal distribution of the particle energy. In the case of the injection and simultaneous 
' escape the resulting spectra depend mostly on the energy distribution of the injected particles. In the simplest case, where the particles are 
t| injected at the lowest possible energies, the competition between the acceleration and the escape forms a power-law energy distribution. 
We apply our modeling to the high energy activity of the blazar Mrk 501 observed in April 1997. Calculating the evolution of the electron 
spectrum self-consistently we can reproduce the observed spectra well with a number of free parameters that is comparable to or less than in 
• • the "classic stationary" one-zone synchrotron self-Compton scenario. 
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1. Introduction 

The overall emission from blazars (e.g. Urry & Padovani l 19951 1 
is successfully explained as the result of synchrotron and 
inverse-Compton emission by relativistic electrons within a 
relativistic jet (e.g. Ghisellini et al. 119981 for a different view 
see e.g. Mucke l2003t . Most of the models usually applied to 
describe the emission adopt a phenomenological view, assum- 
ing that some (usually unspecified) mechanism is able to pro- 
duce the electron distribution that is subsequently injected into 
the emission region (e.g. Mastichiadis & Kirk . 1997. Chiaberge 
& Ghisellini|l999 Kataoka et al. 20011 Moderski et al.^OOS I. 
Two processes are usually considered as a possible source of 
the relativistic particles: the first-order Fermi acceleration at 
a shock front and the second-order Fermi acceleration by a 
plasma turbulence in e.g. a downstream region of the shock. 

In the first-order process, the particle crossing the shock 
discontinuity gains an amount of energy that is proportional 
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to the shock velocity. Behind the shock the particle may be 
scattered back by magnetic inhomogeneities, and as a result 
it crosses the shock many times, each time gaining energy. 
This scenario is well known as a relatively efficient accelera- 
tion process that provides a power-law energy spectrum (e.g. 
Bell ,1978. Bland ford & Ostriker [T979l Drurv [T9"83. Blandford 
& Eichler [79871 Jones & Ellison n'991t . This power-law spec- 
trum is generated by a kind of competition between the accel- 
eration and the escape of the particles from the shock region. In 
a simple approach, the shock acceleration may be described in 
a deterministic way, where the efficiency of the process is de- 
scribed by the characteristic acceleration time (f^co e.g. Kirk et 
al. ll998> . With this approach the particle energy gain per time 
unit is precisely defined. 

The second-order Fermi acceleration assumes reflection of 
the charged particles by a magnetized cloud or, in more so- 
phisticated cases, by the magnetic inhomogeneities or plasma 
waves. The particle may gain or lose energy depending on 
whether the "mirror" is approaching or receding. However, the 
probability of the head-on collisions is higher than for the rear- 
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on reflections. Therefore, on average, the particle can gain en- 
ergy. The main energy gain per bounce is proportional to the 
square of the mirror velocity. However, the net energy gain de- 
pends not only on the mirror velocity but also on the scattering 
rate. Note that the net energy gain for the rear-on collisions is 
larger than the energy losses (e.g. Longair'1992 i. The second- 
order Fermi acceleration is usually treated as a stochastic pro- 
cess and described as the diffusion of the particle energy. This 
process has been successfully applied in many diff'erent astro- 
physical sources (e.g. Eilek & Henriksen 1984 Schlickeiser 
[T9PI & IT989] Dung & Petrosian 1994 Miller & Roberts 1995 
Dermer et al.'1996 Petrosian & Liu 2004 1. However, the effi- 
ciency of the second-order process was quite frequently ques- 
tioned in comparison with the first-order process. Therefore 
the turbulent acceleration was frequently neglected when inves- 
tigating the electron spectrum evolution in extragalactic jets. 

The numerical simulations performed recently by Virtanen 
& Vaino ( 2005 1 show that the efficiency of the second-order 
Fermi acceleration may be comparable to the efficiency of 
the shock acceleration. Therefore, the turbulent acceleration 
may significantly affect the evolution of the particle spectrum. 
Moreover, the acceleration process may be quite complex. For 
example, in the first step, the particle may be efficiently ac- 
celerated at the shock front. After the escape into the down- 
stream region of the shock, the particle is still accelerated by 
turbulent plasma waves. Therefore, the particle energy may in- 
crease enough to let the particle re-enter the shock accelera- 
tion region. In such a case the energy spectrum formed by the 
stochastic process may be re-accelerated by the shock. 

For the sake of simplicity, we do not describe the acceler- 
ation scenario in details. We assume that the acceleration pro- 
cess has a stochastic nature and describe it as the diffusion of 
the particle momentum. This assumption does not mean that 
we consider only the second-order Fermi acceleration. We as- 
sume that the acceleration scenario is quite complex as in the 
example described above. Therefore the particle energy gain 
should be treated as a stochastic process and described by the 
diffusion. The efficiency of the diffusion is a free parameter in 
our model. We analyse the influence of this process in detail for 
the evolution of the particle spectrum. Moreover, we also anal- 
yse the evolution of the emission generated by the accelerated 
electrons inside a homogeneous spherical source. The electrons 
are producing the synchrotron and inverse Compton (hereafter 
IC) radiation. The synchrotron radiation field inside the source 
is scattered by the same population of electrons that produces 
this radiation. This is the synchrotron self-Compton process 
(hereafter SSC). This simple model may have direct applica- 
tion to many astrophysical sources and in particular to those 
BL Lac objects, that produce the very high energetic gamma 
rays observed up to the TeV band (e.g., Catanese et al. fl997l 
Maraschi et al l 1999 ' Krawczynski et al. l200H . 

In all the simulations performed in this work, the accelera- 
tion process is fully compensated for at some energy by the ra- 
diative cooling. This compensation leads to the stationary dis- 
tribution of the particle energy. In a first approach we analyse 
a few evolutionary scenarios, where we neglect the possible 
injection or escape of the particles during the evolution of the 
source. Then we investigate a few cases where we introduce the 



continuous injection and escape of the particles from the accel- 
erating region. Our results can be used to investigate the origin 
of the variability of Mrk 501 observed in April 1997 (Plan et al. 
1998 Djannati-Atai et al. 1999 1. We wifl show that our model 
successfully explains the different high states observed in this 
source, without having to use any more free parameters than 
the most commonly used one-zone SSC model. 

2. The momentum-diffusion equation 

The main difference between this paper and the previous works 
about modeling the TeV blazar emission is in the description 
of the particle acceleration process. We describe it as diffusion 
in the particle momentum space, where the evolution of the 
isotropic, homogeneous phase-space density (/) is described 
by the momentum-diffusion equation: 



dp 



dt 



P Dip, t) 



dp 



(1) 



where p - fSyis the dimensionless particle momentum, D{p, t) 
the momentum-diffusion coefficient, y the particle Lorentz fac- 
tor (p - Vt^^^^)' ^iid P the particle velocity in units of c (e.g. 
Borovsky & Eilek ^86 1. The particle number density (A^) is 
directly related to the phase-space density: 



N{p,t)^A7ip'f{p,t). 

Therefore, we can rewrite this equation: 
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describes the average efficiency of the acceleration process. 
In the case of ultra-relativistic particles (J3 - I), the particle 
momentum becomes equivalent to the particle Lorentz factor 
(p = -y). Therefore we can rewrite the equation in the form: 



dNjy, t) 
dt 
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(e.g. Brunetti f2004> . In order to describe the radiative cooling 
of the particles and their possible escape or injection, we have 
to introduce three more terms 



dNiy, f) 



d_ 

dt dy 

- E(y,t) + Q(y,t) 
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(6) 



The C(y, t) parameter describes the synchrotron and IC cooling 
of the particles 



C(r,o = T 



4 o-jc T r 

r [f/B(0 + f/md(r,0], 



3 nipC^ 



(7) 



where U b and f/,ad are magnetic and radiation field energy den- 
sities, respectively, is the electron rest mass, and the ctj is 
the Thomson cross section. The escape of the particles is de- 
scribed by the E{y, t) term. For the sake of simphcity, in all 
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calculations presented in this work we assume that the escape 
term is independent of energy: 



N(y,t) 



(8) 



characterized by a constant escape time (fesc) related to the 
source size. The injection of the particles is described by Q(y, t) 
and we define this function separately for each test. 

Finally, to complete the description of the equation, we 
have to specify the diffusion coefficient. In all our tests we as- 
sume a Fermi-like acceleration process and define the coeffi- 
cient as a time-independent power-law function of the particle 
energy 



D{j) = |r', 



(9) 



where = 1/facc represents the efficiency of the acceleration 
process described by a characteristic acceleration time (face)- 
The acceleration time is a free parameter in our model. The 
above definition of the diffusion process gives the following 
relation for the acceleration term 



Aiy) = ^. 

tare 



(10) 



Note that this acceleration term describes the particle energy 
gain per time unit precisely and is identical to the terms used 
by other authors (e.g. Kirk et al. I1998> in those simulations 
where the acceleration process is described in a deterministic 
way. This helps to directly compare our test with the other sim- 
ulations. 

3. A quasi-thermal stationary solution 

In the first approach we analyse a stationary solution (A^ = 0) 
of Eq.|6lin the case of the simplified cooling (C = Coy^, i.e. no 
synchrotron self-absorption is considered, and the IC cooling 
is assumed to occur entirely in the Thomson regime) and with- 
out injection and escape of the particles during the evolution of 
the system 



d_ 

dy 



{C{y)-A{y))N{y) + D(y) 



dN(y) 



dy 
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(11) 



Note that the above assumption requires that the initial energy 
distribution N(y, f = 0) = Ni„i{y) + 0. According to Chang 
& Cooper J1970> . the general solution of this equation has the 
form 



f^(y) - jcexp 



■J D(r') 



(12) 



where x is an integration constant. For our particular assump- 
tion about the diffusion process, the stationary solution is an 
ultrarelativistic Maxwellian 



f^(y) — X exp 



^^dy 



xy exp[-2Cofacc(r - 1)1- 



(13) 



The maximum of this function appears at the equilibrium en- 
ergy (7e), where the acceleration process is fully compensated 



for by the cooling. If we describe the cooling using the charac- 
teristic cooling time 



Aiy) = 



1 



Cor' 



then at the equilibrium we have fcooi(Te) - face that gives 



1 



face Co 



(14) 
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This simple solution shows a fundamental difference be- 
tween the evolution where the acceleration is treated as a dif- 
fusion process and the acceleration described only by the de- 
terministic acceleration term. In the case of the deterministic 
acceleration, the system will try to reach the stationary state 
accelerating all the particles to the equilibrium energy and to 
build a monoenergetic population of the particles at this energy 
("pile-up"). 

Let us analyse now a more general case where the diffusion 
and cooling coefficients are defined with power law functions 
{D - DqY' and C = Coy''). The stationary solution for this case 
is given by 



N(y) - xy exp 



Co -l) 



Do(fl -\-b) 



(16) 



where the power law part of the function always has the same 
slope (y^), but the shape of the exponential part depends on the 
relation between the particular properties of the diffusion and 
cooling processes. Therefore in some physical scenarios, we 
expect the formation of a quasi-Maxwellian distribution, not a 
perfect Maxwellian. It should also be mentioned here that in the 
stationary spectrum only the value of the constant x depends on 
the initial energy distribution. If we know the total number of 
the particles in the system 



Mr 



Nmi(y)dy, 



then we can easily find the value of this constant 



y^ exp 



Co(y'»'-°-l) 



Do(a-\-b) 



(17) 



(18) 



dy 



The most important conclusion that appears from the anal- 
ysis of the stationary solution is that the acceleration of the par- 
ticles described as the diffusion process when competing with 
the radiative cooling may lead to a thermal or quasi-thermal 
particle distribution. This effect was showed for the first time 
by Schlickeiser ( 1984i in a somewhat more complex way. The 
temperature of this thermal or quasi-thermal distribution is as- 
sociated to the value of the equilibrium energy. The shape of 
the stationary spectrum is independent of the initial distribu- 
tion, which only determines the absolute level of the stationary 
solution. 

Recently Sauge & Henri (12004) proposed a scenario where 
the quasi-Maxwellian distribution of the electrons obtained 
from the acceleration, described as a diffusion process, is in- 
jected into a spherical region where the particles are producing 
the SSC emission. They successfully apply this model to the 
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activity of Mrk 501 observed in April 1997 (Plan et al. [T^ 
Djannati-Atai et al. ll999l l. However, they did not study the ac- 
celeration process by assuming an arbitrary spectrum of the in- 
jected particles. 

A log-parabolic energy distribution, which is similar to 
a thermal distribution, has been proposed by Massaro et al. 
J2004> in order to explain the X-ray and TeV spectra of Mrk 
421. However, this log-parabolic distribution was obtained by 
assuming some specific nature of the acceleration process. 

4. Numerical approach 

In several simple cases, it is possible to find an analytic, time- 
dependent solution to Eq. |6l However, in the case of the SSC 
emission, where the radiation field is interacting with the par- 
ticles that have generated this radiation, this partial differential 
equation becomes a partial differential integral equation. This 
poses a large problem for a attempt to analytically solve this 
equation. Therefore, in our simulations we use the numerical 
method to find the time-dependent solutions. 

The partial differential equation that describes the evolu- 
tion of the particle energy is a Fokker-Planck type equation. 
We use a very useful numerical differentiation scheme in our 
computations proposed for such type of equations by Chang & 
Cooper (1970 1. This method assumes the forward differentia- 
tion in time and the centred differentiation in the energy that is 
continuously transformed into the forward differentiation (see 
e.g. Press et al. 1989 for more information about the numerical 
differentiation). The transformation is correlated with the cur- 
rent position of the equilibrium between the acceleration and 
cooling. This is the big advantage of this method that provides 
stability for the numerical solution in the presence of equilib- 
rium. Note that a constant differentiation method in the en- 
ergy range (e.g. centred) can be used effectively only in the 
case of a monotonic evolution of the spectrum towards lower 
(face » fcool) or higher (face feooi) energies during the evolu- 
tion of the system (e.g. Chiaberge & Ghisellini 1999 Moderski 
et al. l2003> . Moreover, the constant differentiation method also 
limits the shape of the initial spectrum or the injection func- 
tion. The relatively "sharp" initial/injected function (e.g. Dirac 
delta function) may cause a large numerical diffusion, which 
can be reduced if we increase the number of the mesh points at 
the cost of calculation efficiency. However, this increase may 
quickly lead to large numerical instabilities (see e.g. Press et 
al. .l989.) . The method that we use is free from these problems. 
It is an implicit method, which means that it is necessary to 
solve a system of linear equations in order to obtain the solu- 
tion. This provides additional stability for the numerical com- 
putations. The system of equations forms a tridiagonal matrix 
that can be easily solved numerically (e.g. Press et al. ll989> . In 
our computations we have slightly modified the originally pro- 
posed method by using logarithmic steps in the energy range. 
This helps to greatly reduce the number of mesh points neces- 
sary to obtain a precise result. 

The additional advantage of this differentiation scheme is 
that the method provides conservation of the total particle num- 
ber in the absence of external sources or sinks. Therefore, the 
result of numerical computation can be compared, for exam- 



ple, with the analytic stationary solutions discussed in the pre- 
vious section. For the typical number of the time and energy 
mesh points used in our computations, we obtained conser- 
vation of the total particle number with precision better than 
0.01%. Moreover, we tested our numerical code with some ba- 
sic time-dependent analytical solutions and found very good 
agreement. 

5. A simple evolution of the energy spectrum 

In this section we present a few simple cases the energy spec- 
trum that evolves in the case of the stochastic acceleration. In 
all tests we assumed some initial energy distribution and chose 
a constant efficiency for the acceleration process (A = y/faee) 
that in the competition with a constant cooling (C = Coy^) pro- 
vides equilibrium at an energy equivalent to ye = 10'*. 

5.1. No injection and escape 

In the first three tests we neglected injection and escape of the 
particles during the evolution of the system, and we analysed 
only the evolution of the initial spectrum. 

• First, we assume a very naiTow initial distribution (1 < y < 
2) with an average energy of the particles significantly less 
than the equilibrium energy. The density of the particles is 
constant over the whole energy range of the initial spec- 
trum. In this particular case t^^c ^ feooi for the initial en- 
ergy of the particles. Therefore, the acceleration process 
tries to increase the energy of all the particles. In the case of 
the stochastic acceleration, some of the particles may reach 
the equilibrium energy or even a higher energy after an in- 
finitesimally short time of the system evolution (fo -i- At). 
However, the number of the particles that can gain this en- 
ergy is negligibly small in comparison to the initial num- 
ber of particles. The system needs more time to accumu- 
late a significant number of the particles around the equi- 
librium energy. The simulation presented in Fig.^ shows 
a systematic decrease in the particle density for the initial 
energy range and a simultaneous increase in the density 
around the equilibrium energy. In this particular case the 
densities in these two energy ranges become comparable 
after the evolution time of fgyo ^ 5 face- After a longer time 
(fevo S 20 face), the acceleration will significantly decrease 
the density of the particles with the initial energy, accumu- 
lating a dominant part of the particles around the equilib- 
rium energy and reaching the stationary Maxwellian distri- 
bution. In the stationary state all of the energy acquired by 
the particles in the acceleration process is radiated away. 
Note that the deterministic acceleration in this particular 
test could shift the initial spectrum towards the equilibrium 
energy without significantly modifying the shape as long 
as the average energy of the particles in the system is sig- 
nificantly less than ye- Close to the equilibrium any initial 
spectrum should be transformed by the deterministic ac- 
celeration into a quasi-monoenergetic distribution. Finally, 
the deterministic acceleration appears less efficient than the 
stochastic acceleration. For the same value of the accelera- 
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Fig. 1. An example of the evolution of the electron spectrum 
in the case of no injection and escape of the particles. The up- 
per panel (a) shows the evolution of the initial spectrum (bold 
dashed line) where the particle energy (1 < y < 2) is signifi- 
cantly lower than the equilibrium energy (ye - 10*). The oppo- 
site case, where the initial energy (10^ < j <2x 10'') is higher 
than the equilibrium energy, is presented in the middle panel 
(b). The lower panel (c) shows a case where the lower energy 
in the initial spectrum is lower than ye, but the upper initial en- 
ergy is higher than the equilibrium energy (10 < y < 10^). The 
stationary Maxwellian spectrum on all figures is indicated by 



tion time, this process requires at least twice as much time 
to build a particle density that is comparable to the density 
provided by the stochastic acceleration around the equilib- 
rium energy. 

• In the second test (Fig.^^) we again assume a very narrow 
initial distribution, but now the average energy of the par- 
ticles is significantly higher (10^ < y < 2x 10^) than the 
equilibrium energy. At the beginning the cooling process 
dominates the evolution of the particle spectrum (fcooi ^ 
face). Therefore the initial distribution moves towards the 
equilibrium energy. Note that in the case of a determin- 
istic acceleration this particular initial spectrum could be 
transformed into a quasi-monoenergetic distribution. This 
is related to the fact that the cooling time is shorter for 
the particles with the highest energy. In our test, instead of 
the monoenergetic distribution, we obtain a kind of peaked 
spectrum that is more extended in the energy range than 
the initial spectrum. However, the peak energy of this dis- 
tribution evolves in the same way as the energy of the mo- 
noenergetic distribution in the case of no diffusion. When 
the peak energy reaches ye, the shift in the peak ceases be- 
cause the cooling process becomes fully compensated for 
by the acceleration. On the other hand, after this moment in 
the evolution, the left part of the peaked distribution starts 
to change slope from an exponential to a power law. When 
the index of this power law becomes equal to 2, this evolu- 
tion stops and the formation of the stationary Maxwellian 
distribution is completed. 

• In the next test (Fig. ^), the assumed initial spectrum is 
relatively broad (10 < y < lO*"). The initial minimal energy 
of the particles is significantly less than the equilibrium en- 
ergy, but the maximal energy is significantly higher than ye. 
For the sake of simplicity we also assume a constant parti- 
cle density in this energy range. The specific assumption 
about the initial spectrum makes this case more complex 
than the previous ones. However, the evolution in this case 
in some specific energy ranges is very similar to the cases 
just discussed above. The change in the distribution before 
the minimal initial energy is similar to the evolution pre- 
sented in our second test below the peak energy. Above the 
minimal energy, but below the equilibrium, the changes in 
the spectrum are similar to the evolution presented in our 
first test. Finally, above ye the evolution is analogous to the 
variations presented in our second test, above the equilib- 
rium. 



The tests described above show that indeed the stationary 
state is almost independent of the initial distribution in the ab- 
sence of the injection and escape of the particles. Note that a 
constant escape of the particles {E{y,t) - N{y,t))ltesc) could 
produce a systematic decrease in the particle density propor- 
tional to a factor exp(-f/fesc) without modifying the spectral 
shape (e.g. Schlickeiser ( 1984|. However, the escape and a si- 
multaneous continuous injection are able to completely change 
the evolution and the resulting stationary spectrum. We discuss 
such scenarios in the next subsection. 
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5.2. Continuous injection and escape 

In the next step we discuss three more cases where we assume 
continuous injection and escape of the particles. The escape in 
our calculations may be explained in two different ways. First, 
this term may describe a real escape of the particles into the 
region of a source where the magnetic field strength is signifi- 
cantly smaller and therefore the efficiency of the particle emis- 
sion is also significantly less. The escape may also approximate 
adiabatic losses in the particle energy. The description of the 
adiabatic cooling inside an expanding source is slightly differ- 
ent from the description of the escape. However, the precise 
description of the SSC emission of an expanding source is be- 
yond the scope of this paper. 

In three more tests presented in this subsection we assume 
that the profiles of the continuous injection are equivalent to 
the initial distributions used in the previous tests: 

• The injection profile in our fourth test is identical to the 
initial distribution in the first test (Fig. |2t). The particles 
are continuously injected at the lower energy (1 < y < 2) 
and systematically accelerated up to the equilibrium en- 
ergy. The escape of the particles is described by the charac- 
teristic escaping time fesc = face- The competition between 
the acceleration and the escape in the case of continuous in- 
jection produces a power law distribution that extends from 
the maximal energy of the injected particles up to ye- Below 
the maximal injected energy the spectrum is almost flat. 
Above the equilibrium energy the diffusion process pro- 
duces an exponential cut-off. This makes this distribution 
different from the spectrum that could be produced in this 
case by the deterministic acceleration, where a pile-up or 
a roll-over is produced before the equilibrium energy, de- 
pending on whether the particle slope is softer or harder 
than 2. Moreover, the slope of the spectrum produced by 
the stochastic acceleration 



1 + 



2fe: 



(19) 



differs from the slope that could be obtained from the de- 
terministic acceleration 



n - 1 + 



(20) 



(e.g. Kirk et al. ll998> . Note that basic shock-acceleration 
models postulate that face - fesc, or in more sophisticated 
cases, postulate that f^cc and fesc have the same energy de- 
pendence. Therefore the index of the particle energy spec- 
trum predicted by such models should be « ^ 2. In our 
acceleration scenario the value of the characteristic accel- 
eration time is a free parameter, and the resulting index 
depends on the initial assumption. Moreover, face and fesc 
may have a different energy dependence in a more com- 
plex scenario. This may give much more complex spectra. 
However, as long as the energy dependence of the charac- 
teristic acceleration and escape time is defined by a power- 
law function, the resulting energy spectrum should also be 
a power-law function. 



o 



T 1 r- 
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0.02 








0.02 
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Fig. 2. An example of the evolution of the electron spectrum in 
the case of continuous injection and escape of the particles. The 
bold dashed line shows the initial spectrum and the profile of 
the injection. We use the same parameters for the initial spec- 
trum and the injection as in the simulations presented in Fig.^ 
For a comparison, the stationary Maxwellian distributions ob- 
tained in the previous tests are presented by a thin dashed line. 



The continuous injection above the equilibrium energy pro- 
duces a power law spectrum up to the equilibrium. The 
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evolution in this energy range is dominated by the cool- 
ing and escape process and therefore the index of the pro- 
duced power law is equal to -2. Below the peak energy, the 
spectrum is systematically changed from an exponential to 
a power law at the stationary state. This change is similar 
to the evolution presented in our second test. However, the 
index of the resulting power law 

„:^-2-^, (21) 

this time is steeper than the slope necessary to obtain the 
Maxwellian spectrum in the second test. 
• In the last test the particles are injected below and above 
the equilibrium energy (Fig.|2l;). The evolution below the 
minimal injected energy is analogous to the evolution in 
the previous test below ye- Above the minimal energy, the 
process is dominated by the injection and by the escape that 
leads to a flat stationary spectrum. Above the equilibrium 
energy the evolution is almost the same as in the previous 
test. 

The simulations presented in this section show that the in- 
jection and escape may dominate the evolution of the spectrum 
and may lead to completely different stationary states. Note that 
in some of our simulations we assume for the sake of simplicity 
a very low energy of the initial or injected particles (1 < y < 2). 
However, the turbulent acceleration of such low energetic lep- 
tons may be diflicult in an electron-proton plasma where ther- 
mal protons very efficiently damps resonant Alfven waves with 
the frequency equal to their gyrofrequency. Also a shock wave 
may have difficulty with the acceleration of such low energetic 
leptons. Only the particles with a Larmor radius larger than 
the thickness of the shock are able to 'feel' the shock discon- 
tinuity. Therefore, some pre-acceleration of the particles is re- 
quired in order to let the turbulent or shock acceleration operate 
efficiently. Such pre-acceleration in our test could be simply 
taken into account by assuming higher energy for the initial or 
injected particles (e.g. 10 < y < 20). Such a relatively small 
energy increase has a negligible impact on the evolution of the 
particle spectrum at very high energies (y ^ lO*"). Therefore, 
for the sake of simplicity we neglect the problem of the pre- 
acceleration in the tests, where we assume a relatively low en- 
ergy for the initial or injected particles. 

6. Evolution of the SSC emission 

The evolution of the synchrotron emission depends directly on 
the evolution of the electron spectrum that also controls the 
inverse-Compton scattering process. We analyse here the SSC 
emission generated by a homogeneous, spherical source filled 
by electrons and a uniform magnetic field. The particles are ac- 
celerated uniformly throughout the emitting region. Out of the 
6 evolutionary scenarios described in the previous section we, 
only investigate here two possible evolutions where the initial 
energy of the particles is relatively low (1 < y < 2). Using 
this simple model we try to explain the high-energy activity of 
MrkSOl. 



6.1. No injection and escape 

In the first test we assume no injection and escape of the parti- 
cles from the acceleration region. This simulation is very simi- 
lar to the first evolutionary scenario discussed in Sect.|5l 

The main difference between the present simulation and the 
previous one appears in the particle cooling conditions. In the 
previous test we assumed, for the sake of simplicity, a con- 
stant cooling (C = Coy^), which can be interpreted as the syn- 
chrotron emission in a constant magnetic field (Ub » f/rad)- 
Here, since we apply our results directly to the emission of the 
TeV blazar Mrk 501, we cannot neglect the radiative cooling 
due to the IC scattering. Moreover, we have to take the de- 
crease of the scattering efficiency in the Klein-Nishina regime 
into account. This effect is approximated in our calculations by 
a constraint on the radiation-field energy density used for the 
computation of the cooling 

f/rad(r,0-— J hy.(v,t)dv, (22) 

where I^yn is the intensity of the synchrotron emission, Vx(y) - 
min[vn,ax, 3mj,c^/(4/!y)], and v^in, v^ax are the minimal and 
maximal frequencies of the synchrotron emission, respectively. 
According to this formula for the very high energetic particles 
(y k, 10''), we have Uy^d S in this particular test. On the 
other hand, the radiation field energy density is much larger 
than Ub for the low energy particles (y < 10^). This, as pointed 
out by Dermer & Atoyan 2002 provides a quite complex form 
of the cooling term, which introduces a significant difference 
in the evolution of the electron spectrum (see also Moderski 
et al. 2005 1. This difference is very visible if we compare the 
stationary state in the previous test (Fig.[^) and the stationary 
spectrum in the present simulation (Fig.^t). Previously we ob- 
tained a pure, Maxwellian distribution whereas now we have 
two bumps in the spectrum. The high-energy bump, that dom- 
inates the spectrum appears at the equilibrium energy related 
mostly to the acceleration process and the synchrotron cooling. 
The bump that appears at the lower energies is controlled by 
the acceleration and the IC cooling that is very strong only for 
the low energy particles. However, as long as the amplitude of 
the low energy bump is smaller than the amplitude of the high 
energy bump, the existence of the low energy bump is not vis- 
ible in the synchrotron spectrum (Fig.^')- The simple formula 
(a - (n - l)/2) that describes the relation between the electron 
spectral index (n) and the index of the synchrotron emission (a) 
becomes invalid when n < 1/3. In such a case the spectral in- 
dex of the power law part of the synchrotron emission becomes 
constant and a = -1/3. This value is equivalent to the spectral 
index of the synchrotron emissivity of a monoenergetic popula- 
tion of the electrons, since in this case the electrons belonging 
to the high energy bump dominate the total synchrotron emis- 
sivity. 

The evolution of the synchrotron radiation presented in 
Fig. I^J" shows that the vF{v) peak of the emission is system- 
atically increasing the amplitude and moving towards higher 
frequencies. However, the movement towards higher frequen- 
cies slows down significantly when the maximal particle en- 
ergy becomes comparable with the equilibrium energy. On the 
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Fig. 3. The upper panel (a) shows the evolution of the electron 
spectrum inside a homogeneous, spherical source in the case of 
no injection and escape of the particles. The remaining panels 
show the synchrotron (b) and the inverse Compton (c) radiation 
of the source. In this particular test we applied our model to the 
high energy emission of Mrk 501 observed almost simultane- 
ously by the BeppoSAX (panel-b, Pian et alTWl and CAT 
(panel-c, Djannati-Atai et al. 1999 1 experiments in April 1997. 
The lower level of the emission (7* of April) is reproduced 
well by the source emission at the evolution time equal to 4.5 
/acc while the highest level (16* of April) appears to be repro- 



other hand, the amplitude increases even if the peak frequency 
appears almost constant. The spectral index below the peak fre- 
quency changes simultaneously with the increase in the peak 
amplitude. However, as discussed above, the limiting value for 
the index is a = -1/3, and this slope remains constant even if 
the peak amplitude is still increasing. 

The changes in the IC emission, presented in Fig. 
are very similar to the evolution of the synchrotron radiation. 
However, the evolution of the peak position appears not to be 
as strong as in the case of the synchrotron emission. This is re- 
lated mostly to the fact that the IC scattering is limited by the 
Klein-Nishina effect. Moreover, the very high energy gamma 
rays (E > ITeV) are significantly absorbed due to interaction 
with the intergalactic infrared radiation field. Since we want 
to apply this simulation to the specific case of Mrk 501, the 
spectra shown take into account the absorption suffered by TeV 
photons interacting with the cosmic infrared radiation back- 
ground, and therefore these cases assume a distance between us 
and the source that is equal to the distance of Mrk 501 . The ab- 
sorption produces a cut-off in the spectra above the frequency 
~ 10^^ [Hz] in this particular case. To the calculate the absorp- 
tion, we use the absorption coefficient derived by Kneiske et al. 
(120021 ED04> . The spectral index of the IC emission below the 
peak is also limited by the value a = -1/3. 

As mentioned above, in order to check whether the dis- 
cussed scenario is able to explain the high energy emission 
of TeV blazars, we apply the results of the simulation to the 
well-known high energy activity of Mrk 501 observed by 
BeppoSAX (Pian et al. 1998 1 and CAT (Djannati-Atai et al. 
[19991 in April 1997. What is important for these data is that 
the X-ray /TeV observations were made almost simultaneously. 
However, we have to point out that there is a delay of a few 
days between the higher (16* of April) and the lower (7* of 
April) emission levels and these two emission levels were ob- 
served during different flaring events (Catanese et al. I1997> . 
Moreover, the X-ray (as well as TeV) spectra were gathered 
during a few hours of integration. The observed variability time 
scales of TeV blazars at X-ray/TeV energies are usually very 
short, from several minutes up to a few hours. 

Our main aim was to show that our proposed model is able 
to reproduce such very short variability time scales with rea- 
sonable assumptions about the physical parameters. We apply 
our instantaneous spectra for the observations integrated over 
some period of time to roughly constrain the physical parame- 
ters. The presented data are one of the best quasi-simultaneous 
observational set ever obtained for TeV blazars. We used the 
same evolutionary scenario to explain both presented emission 
levels. This means that we use exactly the same values for the 
physical parameters in both cases. As shown in Fig.lSji/c, our 
modeling reproduces the observed spectra well if we select the 
radiation generated at the evolution time fgvo = 4.5facc for the 
lower emission level and fevo = 5facc for the higher level. The 
characteristic acceleration time assumed in the calculations is 
equal to R/c, where R = 3.2 x 10'^ [cm] is the source radius. 
We used also the Doppler factor <5 = 21, an initial particle den- 
sity Kini - 7 [cm"-'], and the constant magnetic field intensity 
B = 0.05 [G]. 
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6.2. Continuous injection and escape 

The second test presented in this section assumes continuous 
injection and the escape of the particles. As already shown (Fig. 

this scenario leads to a stationary solution that has a power 
law shape. 

The cooling conditions in the present simulation are sim- 
ilar to the conditions described in the previous subsection. 
However, the complexity of the cooling in the Klein-Nishina 
regime now has no impact on the evolution of the low energy 
part of the spectrum. In the low energy range, the evolution of 
the spectrum is dominated by the acceleration and escape pro- 
cesses. 

In the present simulation, the evolution of the peak position 
is very similar to the evolution in the previous test (Fig. ^p). 
The value of the spectral index below the peak also changes 
fast with the increase in the peak amplitude. However, now the 
limiting spectral index a ^ 0.25 is related to the assumption 
face - fesc = R/c made in this particular test. The assump- 
tion, according to Eq. ^] gives the spectral index of the en- 
ergy spectrum n ^ 1 .5 and the above-mentioned index of the 
synchrotron emission. 

If we compare the IC emission in the previous test (Fig.|3j;) 
with the same radiation in the present calculations, we see that 
the evolution in both cases is very similar The scattering in 
the Klein-Nishina regime and the absorption of the TeV emis- 
sion limit the evolution of the peak. The difference between the 
two simulations appears in the slope of the IC spectrum be- 
low the peak that now is n ^ 0.25. The emission level at the 
end of the simulation (stationary particle distribution) in the 
present (as well as in the previous) simulation is much higher 
than the maximal observed level. However, the spectra selected 
from the evolution sequence before the stationary state provide 
a much better explanation for the observations than the spectra 
that correspond to the stationary state. Moreover, the variability 
of the TeV sources, where a fast increase in the observed flux is 
followed by an almost equally fast decrease (quasi-symmetric 
flare profile), may indicate that the source is almost never able 
to reach the stationary state. Therefore, a significantly higher 
level of the emission related to the stationary particle distri- 
bution is rather unlikely. Note also that relatively small mod- 
ification in the high energy "tail" of the particle distribution 
produces very strong changes in the emission level. 

To check the validity of the present simulation we again 
apply our results to the high energy activity of Mrk 501 ob- 
served in April 1997. We can reproduce the observed spec- 
tra well assuming fevo - 4.7facc for the lower emission level 
and fevo - 5.6facc for the higher level. This may suggest that 
in comparison to the previous test, the system needs twice as 
much time to evolve from the lower to the higher emission 
level. However, the source radius assumed in the preset test 
(R = 10^5 cm) is more that three times smaller than the radius 
in the previous simulation {R - 3.2x lO'*" cm). This means that 
the discussed time was even shorter than in the previous case. 
We use the following values of the other model parameters: 
6 = 33, B = 0.1 1 [G], and the the injected particles distribution 
Qinj = 5 X 10"^ [cm"-' s"']. The parameters used in this test 
are quite similar to the parameters used in the previous calcu- 
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Fig. 4. The evolution of the electron spectrum (a) and the SSC 
emission (b, c) in the case of continuous injection and escape of 
the particles. We apply this modeling again for the high energy 
activity of Mrk 501 observed in April 1997 (Plan et al. [T998l 
Djannati-Atai et al. J999J . 



lations. The specific assumptions about some of the parameters 
and the number of free parameters required by the model are 
discussed in the next section. 
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7. Number of free parameters 

In order to explain the SSC emission of a homogeneous spheri- 
cal source filled by a uniform magnetic field and monoenergetic 
population of the electrons, we have to specify the size of the 
source (R), the magnetic field intensity (B), the particle density 
(Ky), the energy (y) of the particles, and the Doppler factor (6) 
if the source velocity is relativistic. This gives only five free 
parameters. However a monoenergetic population of the parti- 
cles generates a synchrotron spectrum with the index a = -1/3 
and the exponential cut-off that do not provide a good fit to the 
observed spectra of TeV blazars. 

Assuming a power-law distribution of the particle energy, 
we have to specify the density (Kmin) of the particles with the 
minimum energy (ymin - 1) instead of Ky and the maximal 
particle energy (ymax) instead of y. Moreover, we also have 
to specify the slope (n) of the electron distribution. This gives 
one more free parameter in comparison to the first model. The 
power-law distribution can explain some of the observations, 
but to explain most of the data we have to apply a more com- 
plex distribution. 

Using a broken power-law energy distribution, we can ex- 
plain most of the observed spectral shapes. However, to de- 
scribe the broken power-law, we have to specify two more 
free parameters. This approach requires additional information 
about the break energy (ybrk) the slope (112) of the sec- 
ond part of the energy spectrum. Therefore the model requires 
eight free parameters. It should be mentioned that in some cases 
when the second part of the electron spectrum is relatively steep 
(«2 » 3), the synchrotron and IC peaks are related to the break 
in the energy spectrum. Therefore y^ax is unimportant, which 
reduces the number of free parameters to seven. However, this 
simplification does not apply to the high-energy observations 
of Mrk 501 used in this work. To explain the data with the bro- 
ken power law and the simple homogeneous model we have 
to specify at least eight free parameters (e.g. Katarzynski et al. 

The model presented in this work uses a more complex 
spectra than either a power or a broken power-law. Our spectra 
are calculated self-consistently according to the assumed phys- 
ical conditions. Therefore, instead of the spectral parameters, 
we have to describe the processes that control the evolution of 
the spectrum (acceleration, radiative cooling, initial spectrum, 
or continuous injection and possible escape). In principle, the 
model requires seven free parameters {6, R, B, Ki„i or Qi„j, t^^c, 
fesc, and fevo)- However, if we assume that the escaping process 
is negligible, then the number of free parameters is reduced to 
six. We made such an assumption in the first simulation that has 
been applied to the activity of Mrk 501 (Sec. 16. it . Moreover, 
we assume in this particular test that the acceleration time is 
equal to the light crossing time of the source (R/c), which re- 
duces the free parameters to five. This is a completely arbitrary 
assumption, but it provides reasonable results. To explain the 
lower of the two observed emission levels, the source requires 
an evolution time equal to 4.5 R/c, in the first simulation. This 
gives about 6.5 hours in the observer's frame. However, to get 
from the lower emission level to the higher level, the source re- 
quires only 0.5 R/c, which gives 0.7 h in the observer's frame. 



This time scales are in good agreement with the variability time 
scales observed in TeV blazars (Maraschi et al. ll999l Fossati et 
al. <2000i Aharonian et al.'2002\. 

In the second test performed to explain the activity of 
Mrk 501, we assumed continuous injection and escape of the 
particles. The characteristic escape time was assumed to be 
equal to the crossing time. This assumption describes on av- 
erage the fastest possible escape of the particles. Moreover, we 
assumed that face is also equal to /?/c as in the previous test. In 
fact the acceleration time cannot be longer than R/c because 
the particles could escape from the source without a signifi- 
cant energy gain. On the other hand, the acceleration time can 
be shorter than R/c, which could coiTespond to a more effi- 
cient acceleration. However, face ^ fesc do not provide the re- 
quired slope of the electron spectrum (n ~ 1 .5 ... 2 a ~ 
0.25 . . .0.5). Assuming face = fesc - R/c we reduce the num- 
ber of free parameters to five (6, R, B, 2inj, fevo)- This gives 
a significant advantage in comparison to a simple homoge- 
neous SSC model with the broken power law electron spec- 
trum. Moreover, our model describes a time-dependent evolu- 
tion of the source emission. This gives the possibility of elim- 
inating one more free parameter. Suppose that we have obser- 
vations of the source in a stationary state. Then the evolution 
time (fevo) can be eliminated. This reduces the number of free 
parameters to only four. 

The reduction of the number of free parameters with re- 
spect to the commonly used SSC models may have very im- 
portant implications. Suppose that we can describe the emis- 
sion with only three free parameters. Then the values of these 
parameters could be directly derived from the relationship for 
the synchrotron peak frequency, the level of the synchrotron 
emission at the peak and the observed variability time scales. 
This means that we could predict the IC emission from the 
observed X-ray radiation and verify for example the absorp- 
tion of the very high energy gamma-rays by the infrared inter- 
galactic radiation field. However, we again have to stress that in 
constructing our model we have made a few completely arbi- 
trary assumptions (e.g. face - R/ s). These assumptions provide 
quite reasonable results but this is not a proof that our model 
correctly describes the source evolution. Our assumptions re- 
quire independent observational confirmation. Our main aim 
was only to show how some of the free parameters could be 
eliminated from the model. If we neglect the arbitrary assump- 
tions, our time-dependent model requires seven independent 
parameters. This is the number of free parameters also required 
by simple stationary SSC scenario. 

8. Summary 

We have proposed a simple time-dependent model for the SSC 
emission of a homogeneous source, which has direct applica- 
tion for the high energy emission of TeV blazars. The most im- 
portant assumption in our modeling is that the acceleration of 
the particles has a stochastic nature. To explain such a case, we 
use the momentum diffusion equation that describes the evolu- 
tion of the particles energy spectrum. For the sake of simplicity 
we have not described the possible acceleration processes in 
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details but assumed only a constant average efficiency of the 
acceleration that is a free parameter in our model. 

We first analysed the stationary analytic solution of the dif- 
fusion equation. We show that for constant cooling conditions 
and in the absence of particle injection or escape, the competi- 
tion between the stochastic acceleration and the radiative cool- 
ing may lead to a thermal or quasi thermal distribution of the 
particles. This result is in a good agreement with the results of 
SchHckeiser ( iT^S^b . 

In the next step we used numerical solutions of the diffusion 
equation to analyse the time-dependent evolution of the elec- 
tron energy spectrum. In the case of no injection and escape of 
the particles, the resulting stationary spectrum is almost inde- 
pendent of the initial distribution. Then, assuming a continu- 
ous injection and simultaneous escape of the particles, we have 
showed that the stochastic acceleration may lead to different 
energy distributions that depend on the assumptions about the 
spectrum of the injected particles. In the most realistic case, 
when particles are injected at the lowest possible energies, the 
stochastic acceleration in competition with the escaping gener- 
ates a power-law particle distribution. 

Finally, we applied our model for the high energy activity 
of Mrk 501 observed in April 1997. We have used two opposite 
scenarios for the evolution of the electron spectrum. First we 
investigated the case where the escape and injection of the par- 
ticles is neglected. In the second test we instead assumed con- 
tinuous injection and escape. In both cases we have obtained 
satisfactory fits for the observed spectra. We have to point out 
that our model describes time-dependent evolution of the emis- 
sion using the number of free parameters that is comparable to 
or even less than the number required by the very simple SSC 
scenario. 

We point out that a general prediction of this model is that 
in the phase of ongoing preequilibrium stochastic acceleration 
of a population of low-energy electrons (see Fig. l3l4l a.b'). the 
emissivity level of the peak of the synchrotron emission and the 
energy of the peak are correlated, a behaviour already observed 
a few times in blazars (e.g. Fossati et al. 2000, Tanihata et al. 
12004 !. We are going to investigate this problem in detail in our 
further work. We will also investigate some specific relations 
between the X-ray and the TeV light curves that are predicted 
by the model. 

Finally, a straightforward extension of our study is the ap- 
plication to the case of powerful blazars, in which the dominant 
role in the cooling of the high-energy electrons is played by the 
external radiation field of the host quasar In these sources, as 
opposed to the case of TeV BL Lacs, the low-energy branch 
of the IC component is well known from X-ray observations, 
providing a supplementary constraint to the model. 
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